Long-range magnetic order and spin-lattice coupling in the delafossite CuFe02 
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The electronic and magnetic properties of the delafossite CuFe02 are investigated by means of 
electronic structure calculations. They are performed using density functional theory in the gener- 
alized gradient approximation as well Eis the new full-potential augmented spherical wave method. 
The calculations reveal three different spin states at the iron sites. Taking into account the cor- 
rect crystal structure, we find long-range antiferromagnetic ordering in agreement with experiment. 
Contrasting previous work, our calculations show that non-local exchange interactions lead to a 
semiconducting ground state. 
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Long-range magnetic ordering on triangular lattices 
with antiferromagnetic exchange interactions is at the 
focus of continuing interest. This is due to the strong ge- 
ometric frustration experienced by such systems, which 
may lead to a large variety of magnetic states including 
incommensurate and non-coUinear spin arrangements. 
As a consequence, there is usually a complex response to 
magnetic fields, which gives rise, e.g., to magnetization 
steps [H, [1: 9 ■ While this behavior has been observed in 
the trigonal chain cobaltates, where it leads to a striking 
spin dynamics 0, the genuine situation is that of well 
separated triangular layers as they are found, e.g., in the 
delafossite-type compounds ABO2. 

In general, this broad class of materials has aroused 
much interest due to a broad range of exciting physical 
properties [3], including strongly anisotropic very high 
electrical conductivities as, e.g., in PdCo02 and semicon- 
ducting behavior in antiferromagnetic delafossites. High 
optical band gaps in Cu- and Ag-based materials allow 
for simultaneous transparency and p-type conductivity 
and, hence, the development of transparent optoelec- 
tronic devices. This large variety is caused by the stack- 
ing of monoatomic triangular layers within the rhom- 
bohedral structure [J. Edge-sharing distorted oxygen 
octahedra surrounding the B-atoms form BO2 sandwich 
layers, which are linked to the A-atom layers via linear 
0-A-O bonds (for the crystal structure see also Fig. 1 
of Ref. 0). Generically, the A and B atoms are mono- 
and trivalent, respectively. Depending on the chemical 
composition, this opens a zoo of behaviors: for instance, 
if the A+ ion is in a (P configuration good metallic con- 
ductivity is observed as in the case of PdCo02, while if 
it is in a d^'^ configuration, the degrees of freedom domi- 
nating the low-energy physics are due to the B atoms as, 
e.g., in CuCr02 and AgCo02. 

As for many other magnetic delafossite compounds, 
the exact magnetic structure of CuFe02 has long been 
a matter of dispute 0, S|. Using neutron diffraction, 
Mekata et al. were able to distinguish two different mag- 



netic phases below Tni = 16K and Tn2 = UK. They 
are connected with monoclinic and orthorhombic mag- 
netic supercells, respectively, of the undistorted rhom- 
bohedral unit cell with commensurate and incommensu- 
rate collinear arran gem ents of the localized 4.4^b Fe'^+ 
moments 0, [13, 
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Il2| . Observation of a noncoUinear- 
incommensurate phase in magnetic field was taken as 
indicative of possible multiferroic behavior [l^ . which 
was indeed observed in Al-dopcd CuFe02 [l^ . Quite re- 
cently. X-ray and neutron diffraction measurements by 
Ye et al. contrasted the previous observations by re- 
vealing structural distortions accompanying the magnetic 
phase transitions, which eventually lead to a monoclinic 
structure at 4 K jl4j . 

Only few electronic structure calculations for magnetic 
delafossite compounds have been reported in the litera- 
ture [3, Il7|,|l8|. From LDA calculations, Galakhov et 
al. obtained a ferromagnetic state for the rhombohedral 
i?3m structure with a magnetic moment at the Fc site of 
about 0.9 /is, much lower than the experimental value 
[H]. The Fe 3d t2g states were found above the Cu 3c? 
states just at Ep, in disagreement with both photoemis- 
sion data and the fact that CuFe02 is a semiconductor 
with an optical band gap of about 1.15 eV. In contrast, 
LDA-I-U calculations led to a band gap of 2eV and a 
magnetic moment of 3.76 /is . However, the occupied Fe 
3d states were located at about 9eV below the valence 
band maximum and thus much to low [l5l |. Recent cal- 
culations by Ong et al. using the generalized gradient 
approximation (GGA) [2^ resulted in a high-spin state 
with a magnetic moment of 3.78 per Fe and the Fe 3d 
t2g spin-up states below the Cu 3d bands in agreement 
with photoemission and X-ray emission data [l7[. How- 
ever, again a finite optical band gap was arrived at only 
after taking into account electronic correlations within 
the LDA-I-U scheme [17} . 

In the present work we apply the new full-potential 
augmented spherical wave method to study the electronic 
properties of CuFe02. In doing so, we take for the first 
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time the experimentally observed low-temperature struc- 
ture into account. Our calculations resolve open issues 
by revealing i) an antiferromagnetic ground state for the 
monoclinic structure in perfect agreement with the ex- 
perimental data, ii) the opening of a fundamental band 
gap already at the GGA level, i.e. without taking lo- 
cal correlations into account, and iii) the quite unusual 
existence of three different magnetic states of assumed 
ferromagnetic CuFe02, which so far has been seen only 
for elemental iron [l3| . 

The calculations are based on density-functional the- 
ory and the generalized gradient approximation (GGA) 
[20| with the local-density approximation parameterized 
according to Perdew and Wang [2]|. They were per- 



formed using the scalar-relativistic implementation of 
the au gme nted spherical wave (ASW) method (sec Rcfs. 
12, HI,]!! and references therein). In the ASW method, 
the wave function is expanded in atom-centered aug- 
mented spherical waves, which are Hankcl functions and 
numerical solutions of Schrodinger's equation, respec- 
tively, outside and inside the so-called augmentation 
spheres. In order to optimize the basis set, additional 
augmented spherical waves were placed at carefully se- 
lected interstitial sites. The choice of these sites as well 
as the augmentation radii were automatically determined 
using the sphere-geometry optimization algorithm [25| . 
Self-consistency was achieved by a highly efficient algo- 
rithm for convergence acceleration [26|. The Brillouin 
zone integrations were performed using the linear tetra- 
hedron method with up to 1156 and 3180 k-points within 
the irreducible wedge of the rhombohedral and mono- 
clinic Brillouin zone, respectively [24, 27|. For the mon- 
oclinic magnetic supercell, up to 100 k-points were used. 

In the present work, a new full-potential version of 
the ASW method was employed 28|. In this version, 
the electron density and related quantities are given 
by spherical-harmonics expansions inside the muffin-tin 
spheres. In the remaining interstitial region, a repre- 
sentation in terms of atom-centered Hankel functions is 
used [2^. However, in contrast to previous related im- 
plementations, we here get away without needing a so- 
called multiple- K basis set, which fact allows to investi- 
gate rather large systems with a minimal effort. 

The calculations used the crystal structure data by Ye 
et al. [14| . As a starting point, spin-degenerate calcu- 



lations for the rhombohedral structure were performed. 
The resulting partial densities of states (DOS) are shown 
in Fig. [TJ While the lower part of the spectrum is dom- 
inated by O 2p states, the transition metal d states lead 
to rather sharp peaks in the interval from —3 to -1-2 eV. 
In particular, the t2g and Cg manifolds of the Fe 3d states 
as resulting from the octahedral coordination are recog- 
nized. This representation of the partial DOS used a 
local rotated coordinate system with the Cartesian axes 
pointing towards the oxygen atoms, cr-typc overlap of 
the O 2p states with the Fe 3d e„ orbitals leads to the 
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FIG. 1: (Color online) Partial densities of states (DOS) of 
rhombohedral CuFe02. Selection of Fe 3d orbitals in this and 
the subsequent figures is relative to the local rotated reference 
frame, see text. 



contribution of the latter near — 5eV. In contrast, due to 
the much weaker 7r-type overlap of the O 2p states with 
the t2g orbitals, these states give rise to sharp peaks in 
the interval from — O.SeV to just above the Fermi energy. 
The latter falls right into the upper part of the t2g mani- 
fold and Fe turns out to be in a d^ state. In contrast, the 
Cu 3d states are essentially limited to the interval from 
—3 to — leV and thus Cu can be assigned a monovalent 
d^^ configuration in close analogy with the experimental 
findings. In passing, we mention the finite dispersion of 
the electronic bands parallel to F-A, which points to a 
considerable three-dimensionality arising from the cou- 
pling between the layers as has been observed also in 
other delafossite materials 0. 

Since the perfect triangular lattice of the rhombohe- 
dral structure does not allow for long-range antiferromag- 
netic order, subsequent spin-polarized calculations were 
performed for an assumed ferromagnetic state in a spirit 
similar to the previous work by Galakhov et al. as well 
as by Ong et al. [l^, 17 1. From our calculations, three 



different configurations were obtained corresponding to 
a low-spin, intermediate-spin, and high-spin moment lo- 
cated at the Fe site. The respective partial densities of 
states are displayed in Fig. O The total energies as com- 
pared to the spin-degenerate configuration and the local 
magnetic moments are summarized in Tab.]!] In general, 
the observation of three different spin states for magnetic 
ions is very unusual and so far has been reported only for 
elemental fee Fe [l^. According to the partial DOS, the 
magnetic moments of the low-spin and intermediate-spin 
states are almost exclusively carried by the Fe 3d t2g 
states, which show a spin splitting of about 1 and 2cV, 
respectively. In the high-spin configuration this splitting 
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FIG. 2: (Color online) Partial densities of states (DOS) of 
rhombohedral ferromagnetic (a) low-spin, (b) intermediate- 
spin, and (c) high-spin CuFe02. 



TABLE L Total energies (in mRyd per formula unit) and 
magnetic moments (in ^s) for different crystal structures and 
magnetic orderings of CuFe02. 



structure 


magn. order 


A TP 




TflQ 


rhomb. 


spin-deg. 


0.0 






rhomb. 


ferro (LS) 


-16.7 


1.03 


-0.02 


rhomb. 


ferro (IS) 


-12.0 


2.02 


-0.02 


rhomb. 


ferro (HS) 


-19.2 


3.73 


0.21 


monoclinic 


spin-deg. 


-6.0 






monoclinic 


ferro (LS) 


-21.5 


1.04 


-0.02 


monoclinic 


ferro (IS) 


-19.0 


2.08 


-0.02 


monoclinic 


ferro (HS) 


-32.0 


3.62 


0.19 


monoclinic 


antiferro 


-46.0 


±3.72 


±0.08 



increases to « 3.5 eV and the magnetic moment is car- 
ried by both the t2g and Cg states. In addition, due to the 
strong cr-type overlap with the latter, the O 2p states also 
experience a substantial polarization. For the same rea- 
son, the eg moments start to form already in the energy 
interval of the O 2p states leading to distinctly different 
spin-up and spin-down Cg partial DOS. A similar behav- 
ior has been also observed in other high-spin systems 
and termed as the formation of local extended magnetic 
moments [s^l. As is obvious from Tab. [H all three fer- 
romagnetic configurations have energies lower than the 
spin-degenerate situation. However, the high-spin state 
is most stable for the rhombohedral lattice. In summary, 
our calculations not only reproduce both the low-spin and 
high-spin results obtained by Galakhov et al. as well as by 
Ong et al. and explain the differences between their find- 
ings but additionally prove the existence of yet another, 
intermediate-spin state. Furthermore, the high-spin par- 
tial DOS compare very well with the photoemission and 
X-ray emission data [15[ . 

In a second step, the monoclinic structure observed by 
Yc et al. was considered [31 • Note that in this struc- 
ture there is still only one Fe atom per unit cell. Both 
the spin-degenerate and the spin-polarized ferromagnetic 
calculations led to essentially the same partial DOS as for 
the rhombohedral structure. In particular, again three 
different magnetic configurations were found with the lo- 
cal magnetic moments as listed in Tab. U being almost 
identical to those obtained for the rhombohedral struc- 
ture. However, the total energies, also given in Tab. HI are 
generally lower by several mRyd with the largest energy 
lowering occurring for the high-spin state. 

Finally, calculations for the eightfold magnetic super- 
cell proposed by Ye et al. [3] were performed. The 
resulting partial DOS are displayed in Fig. [3] and the lo- 
cal magnetic moments and total energy included in Tab. 
m According to these results, the antiferromagnetic state 
has the lowest energy as compared to all other configura- 
tions. In addition, Fe is found to be in a high-spin state in 
agreement with the neutron diffraction data by Mekata et 
al. @,[3l. Remarkably, a band gap of 0.05 cV is obtained. 
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FIG. 3: (Color online) Partial densities of states (DOS) of 
monoclinic antiferromagnetic high-spin CuFe02. 

Thus the non-local exchange interaction included in the 
GGA leads to the semiconducting ground state once the 
monoclinic structure is correctly accounted for. How- 
ever, the band gap is too small reflecting the well known 
shortcomings of the GGA. Additional inclusion of elec- 
tronic correlations, for instance via the LDA+U method, 
is needed to achieve quantitative agreement with experi- 
ment. 

In summary, our calculations for CuFc02 demonstrate 
that i) taking the monoclinic structure into account re- 
sults in an antiferromagnetic ground state in perfect 
agreement with the experimental situation, ii) a funda- 
mental band gap is opened already within the GGA, and 
iii) there is a quite unusual competition among several 
magnetic states, including three different magnetic states 
of assumed ferromagnetic CuFe02, which so far seems 
not to have been obtained for iron compounds. Con- 
cerning the latter point it is remarkable that the trigo- 
nal environment of Fe'^''" renders its magnetic states very 
close in energy, while under most circumstances Fe'^"'" is 
in the high spin state. Even though the effect is less 
pronounced than in the calcium cobaltates, where the 
environment has a dramatic effect on the spin configura- 
tion of the Co'^"'" ions sO, Slj, this opens a route to the 
observation of spin state transitions in Fe'^^ ions as well. 
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